in 
o 
o 

(N 

5h . 
Ph. 

< 

(N 
(N 



On the Stability and Structural Dynamics of Metal Nanowires 

Jerome Biirki, Charles A. Stafford 

Department of Physics, University of Arizona, 1118 E. 4th Street, Tucson, AZ 85721 



a 
a 

c 

o 
o 



> 

ON 

in 

o 
in 
o 

-I— > 
c3 



i 

C 

O 
o 



X 



April 22, 2005 



Abstract This article presents a brief review of the 
nanoscale free-electron model, which provides a con- 
tinuum description of metal nanostructures. It is ar- 
gued that surface and quantum-size effects are the two 
dominant factors in the energetics of metal nanowires, 
and that much of the phenomenology of nanowire 
stability and structural dynamics can be understood 
based on the interplay of these two competing fac- 
tors. A linear stability analysis reveals that metal 
nanocylinders with certain magic conductance values 
G = 1,3,6,12,17,23,34,42,51,67,78,96,... times the 
conductance quantum are exceptionally stable. A non- 
linear dynamical simulation of nanowire structural evo- 
lution reveals a universal equilibrium shape consisting 
of a magic cylinder suspended between unduloidal con- 
tacts. The lifetimes of these metastable structures are 
also computed. 



1 Introduction 

A macroscopic analysis of the mechanical properties of 
thin metal wires suggests that it might be difficult to 
fabricate wires thinner than a few thousand atoms in 
cross section: Consider a cylindrical wire of radius R 
and length L. The maximum stress that the wire can 
sustain before the onset of plastic flow is ay, the yield 
strength. On the other hand, the surface-induced stress 
in a thin wire is a s /R, where a s is the surface tension. 
If a s /R > ay, one would expect the wire to undergo 
plastic flow and, if L > 2nR, to break up under surface 
tension, as in the Rayleigh instability of a column of fluid 
PP. This estimate gives a minimum radius for solidity, 
Rmin — o~ s / o~y ■ The parameters for several simple metals 
are given in Table Plateau realized as early as 1873 
that this surface-tension driven instability of a cylinder is 
unavoidable if cohesion is due solely to classical pairwise 
interactions between atoms 



A great deal of experimental evidence has accumu- 
lated over the past decade, however, indicating that 
metal wires considerably thinner than the above es- 
timate can be fabricated by a number of different 

Even wires with lengths significantly exceeding their cir- 
cumference were found to be remarkably stable jHHHEB 
H7IH8U19j . indicating that some new mechanism must in- 
tervene to prevent their breakup. 

A clue to the resolution of this problem was pro- 
vided by the observation of electron-shell structure in 
conductance histograms of alkali metal point contacts 
[11111411151116) . Like the surface tension, quantum-size ef- 
fects arising from the confinement of the conduction elec- 
trons within the cross-section of the wire become increas- 
ingly important as the wire is scaled down to atomic 
dimensions. In fact, a linear stability analysis |21II22| 
of ultrathin metal wires within the free-electron model 
found that the Rayleigh instability can be completely 
suppressed in the vicinity of certain magic radii. 

In this article, we argue that surface and quantum- 
size effects are the two dominant factors in the energetics 
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Cu 


210 


1.5 


190 


7.1 


2300 


Ag 


140 


1.0 


154 


7.4 


1900 


Au 


100 


1.3 


257 


13 


5600 


Li 


15 


0.44 


99 


29 


26000 


Na 


10 


0.22 


39 


22 


10000 



Table 1 The yield strength ay surface tension a 3 and 
curvature energy 7 S |5| of various monovalent metals. For a 
wire of radius R < a s /ay, the stress due to surface tension 
exceeds ay, signalling a breakdown of macroscopic elasticity 
theory. The electrical conductance G m i n of a ballistic wire of 
radius R m in = Cs/uy is shown in the rightmost column, in 
units of the conductance quantum Go = 2e 2 /h. Note that 
G/Go is approximately equal to the number of atoms that 
fit within the cross section for monovalent metals. 
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of metal nanowires, that is, metal wires with R < i? m i n - 
We show that much of the phenomenology of nanowire 
stability and structural dynamics can be understood 
based on the interplay of these two competing factors. 

This article is organized as follows: In Sec. El we 
describe our continuum structural model for metal 
nanowires. A linear stability analysis of metal nanowires 
is presented in Sec. El Section El describes the struc- 
tural evolution of a metal nanowire from a random 
initial configuration to a universal equilibrium shape. 
The thermally-activated decay of metal nanowires is dis- 
cussed in Sec. El Some concluding remarks are given in 
Sec. El 

2 The Nanoscale Free-Electron Model 

Guided by the importance of conduction electrons in the 
cohesion of metals, and by the success of the jellium 
model in describing metal clusters j2H], the nanoscale 
free-electron model (NFEM) [2U replaces the metal ions 
by a uniform, positively charged background that pro- 
vides a confinement potential for the electrons. The elec- 
tron motion is free along the wire, and confined in the 
transverse directions. Due to the excellent screening [23 
EHJ in metal wires with G > Gq, electron-electron in- 
teractions can in most cases be neglected. The surface 
properties of various metals can be fit by using appro- 
priate surface boundary conditions |27U28| . 

The NFEM is especially suitable for alkali metals, 
but is also adequate to describe shell effects due to the 
conduction-band s-electrons in other monovalent metals, 
such as gold. The experimental observation of a crossover 
from atomic-shell to electron-shell effects with decreas- 
ing radius in both metal clusters [211 and nanowires ^3 
rtl)| justifies a posteriori the use of the NFEM in the later 
regime. 

A nanowire connecting two macroscopic electrodes 
is an open quantum system, for which the Schrodinger 
equation is most naturally formulated as a scattering 
problem. Transport properties can be obtained from the 
scattering matrix using Landauer-type formulas |24U3UI 
I31| . while cohesive properties require the computation 
of the grand canonical potential of the electrons. The 
later can also be expressed in terms of the scattering 
matrix (21], or calculated semiclassically in terms of 
geometrical quantities and a sum over classical periodic 
orbits, as presented in Sect. 12. ll 

Motivated by the argument presented in Tabled the 
ionic degrees of freedom in the wire are modeled as an 
incompressible, irrotational fluid 22,33 . In the Born- 
Oppenheimer approximation, the electronic free energy 
serves as the potential energy for the ions. The ionic dy- 
namics may then be modeled via a surface self-diffusion 
equation |33| . as presented in Sec. l2.3.Tl or. taking ther- 
mal fluctuations into account, via a classical Ginzburg- 
Landau stochastic field theory as presented in Sec. 
l2~T2l 



2.1 Electronic Energy Functional 

Restricting ourselves to axisymmetric structures, the 
grand canonical potential for the electrons /2 e becomes a 
functional of the radius R(z) of the wire. Using the Weyl 
expansion [35j , Q e can be expressed in terms of geomet- 
rical quantities such as the volume V, surface area S, 
and integrated mean curvature C of the wire's surface, 
plus an electron-shell correction, 

n e [R(z),T] =-uV + o- B S-~/ B C+ [ dzV sheU , (1) 

Jo 

where — lo is the bulk value per unit volume, o~ s is the 
surface tension, 7 S is a curvature-energy density, and 
V s heii(R{z),T) is a mesoscopic electron-shell potential, 
shown in Fig.^ The parameters a s and 7 S , tabulated for 
various metals in Table depend on the details of the 
interaction-dependent surface confinement potential |27l 
I28II32| . but can be taken as phenomenological material- 
dependent parameters (along with u>) in our model. The 
leading-order electron-shell correction is, however, in- 
dependent of the Coulomb interaction ,25 ,26 32], and 
therefore insensitive to the details of the confinement 
potential. 

The geometrical quantities S — J Q L dz dS and C = 
f^dz dC are given by 

dS[R(z)] = 2wR(z)^l + (d z R) 2 , (2) 

and 

where d z = d/dz. 

Approximating the confining potential by a hard wall 
at the surface of the wire, the electron-shell potential 
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Fig. 1 Electron-shell potential V s heii (R, T) at zero and two 
finite temperatures, which correspond respectively to 1000K 
and 2500K for Na. The electrical conductance values of the 
magic cylindrical wires are indicated on the upper axis. 
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V s heii can be expressed in terms of a Gutzwillcr-type 
trace formula 



2.3 Structural Dynamics 



to — 1 v=2u; 



(-^) 



fvw COS p^j 
,,2 r 



(4) 



where the sum includes all classical periodic orbits (v, w) 
in a disk billiard |35|. characterized by their number of 
vertices v and winding number w, L vw — 2vR$,vn(Ttw/v) 
is the length of orbit (v,w), and 6 VW — kpL vw — 
3vtt/2. The factor f vw = 1 for v = 2w, 2 otherwise, 
accounts for the invariance under time-reversal symme- 
try of some orbits, and a vw (T) = t vw / smhr vw (t vw = 
nkFL vw T/2Tp) is a temperature-dependent damping 
factor. 

Vsheii(R) exhibits deep minima as a function of R 
(see Fig. GJ, suggesting that some radii are strongly fa- 
vored, which is confirmed by the stability analysis of 
Sec. Note that room temperature is small compared 
to the Fermi temperature Tp — e^/fc^, (in particular, 
T/T F = 0.008 at T = 300K for Na), so that the finite- 
temperature electron-shell potential is essentially indis- 
tinguishable from its zero-temperature limit at experi- 
mental temperatures. 



2.2 Ionic Energetics 

In the Born-Oppcnhcimcr approximation, the electronic 
energy acts as a potential energy for the ionic back- 
ground. The wire can exchange atoms with the macro- 
scopic contacts via surface self-diffusion, so the grand 
canonical ensemble has to be used for the ionic back- 
ground as well, leading to an ionic grand canonical po- 
tential 

fl a = fl e - ll a Nai (5) 

where Af a = V /V a is the number of positive ions in the 
wire (V a = 'STT 2 /k F is the volume of an atom), and /i a 
is the chemical potential for a surface atom in the wire. 
Using Eqs. the ionic free energy ijSJ) becomes 



fin 



dz 



2wa s R(z,t)^l + {d z R) 2 



TTJ S + VsheH (R, T) 



{u + Ha/V a )V, (6) 



where only the leading-order term in the curvature en- 
ergy is included. The chemical potential fx a is obtained 
by calculating the change in the energy (JJJ with the ad- 
dition of an atom at point zq, fi a (zo) — fl e [R(z) + c5(z — 
z ), T] — fl e [R( z ): T~\ , where c = V a /2-KR(z) is chosen so 
that the volume of an atom is added: 



Ha{z) = -uV a + 



V a f 2a s dC[R{z)} dV she 



2ttR 



dR 



(7) 



2.3.1 Surface self- diffusion Since a large fraction of the 
atoms in a nanowire are on the surface, surface self- 
diffusion is the dominant mechanism of ionic motion 33 . 
The dynamics derive from ionic mass conservation: 

7T dR 2 (z,t) d 



V a dt 



+ — [2icR(z,t)J z (z,t)]=0, (8) 



where the z-component of the surface current density is 
given by Fick's law: 



PsD s 



1 



k B T ^/l + (d z R) 2 dz 



(9) 



Here ps and D§ are the surface density of ions and the 
surface self-diffusion coefficient, respectively. The pre- 
cise value of Ds for most metals is not known, but it 
can be removed from the evolution equation by rescaling 
time to the dimensionless variable r = (psDsTp /T)t. 
For comparison to experimental time scales, one can 
estimate that for quasi-one-dimensional diffusion D s w 
vpsa 2 exp(— Es/ksT), where vp> is the Debye frequency, a 
is the lattice spacing, and E s is an activation energy com- 
parable to the energy of a single bond in the solid. Our 
non-linear dynamical model, Eqs. JJHHJ, differs from pre- 
vious studies of axisymmetric surface self-diffusion |3til 
I37U38) by the inclusion of electron-shell effects [last term 
of Eq. (jTJ], which fundamentally alter the dynamics. 

2.3.2 Thermal fluctuations The diffusive dynamics of 
the previous subsection describe relaxation toward struc- 
tures of lower free energy. Once an equilibrium configu- 
ration (i.e., a local minimum of the free energy) is at- 
tained, however, fluctuations about this configuration 
will dominate the dynamics, limiting the dwell time of 
the system in this local minimum. As shown in Sec. 21 
the equilibrium configurations consist of stable cylindri- 
cal nanowires in diffusive equilibrium with unduloid-likc 
contacts 33 3'<J . We therefore study fluctuations of the 
form 

R{z,t)=Ro + (t>(z,t), (10) 

where Rq is the radius of a stable cylinder of length L. 

The energy © can be expanded as a series in <f>. 
For the magic cylinders, corresponding to minima of 
Vsheii(Ro) (cf. Fig. (IJ, the chemical potential for the 
exchange of atoms between the wire and the contacts 
reduces to 

Ma _ a s 

Keeping only the leading-order terms in d z <p, one gets 
fl a = fl a (Ro) +Ti[4>], where fl a (Ro) is the energy of an 
unperturbed cylinder of radius Rq and 



(11) 



H[<j>] = 



-(d z <j>) 2 + V(<j>) 



(12) 
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Here k = 27T(T s i?o and 

V{4>) = V sheU (R + 0) - V sheU (R ) - ^0 2 . (13) 

The problem of stability of nanowires against ther- 
mal fluctuations can be studied as a one-dimensional 
Ginzburg-Landau scalar field theory, perturbed by weak 
spatiotemporal noise, in a domain of finite extent (see 
|3"4"| and references therein): The fluctuations of the 
nanowire radius 4> are treated as a classical field on a 
one-dimensional spatial domain [0, L\. Its dynamics are 
governed by the stochastic Ginzburg-Landau equation 



(14) 



where £(z, i) is unit-strength spatiotemporal white noise, 
satisfying (^(z 1 ,t 1 )£_(z 2 ,t 2 )) = 5(zi~ z 2 )S(t 1 -t 2 ). In Eq. 
I|I4|I , time is measured in units of a microscopic timescale 
describing the short- wavelength cutoff of the surface dy- 
namics |22j . which is given to within a factor of order 
unity by the inverse Debye frequency . The zero- 
noise dynamics is "gradient," that is, at zero temper- 
ature 4> = SH/d(f>, where H[4>] is given by Eq. (H2J. Eq. 
HI4J) represents a considerable simplification compared 
to the volume-conserving dynamics of Eq. JHJ (which in- 
volves derivatives up to d^R), and makes possible an 
analytical treatment of thermal fluctuations. 



3 Linear Stability of Cylinders 

The linear stability of a structure is determined by 
studying the change of energy induced by a small per- 
turbation: If any one perturbation decreases the energy, 
the structure is unstable, while it is stable if all pertur- 
bations increase the energy. 

The most general perturbation of a cylinder of radius 
i?o and length L is 



R(z,(f>) = Ro+X^2J2 b ^(l) e 



i(qz-\-m4 



(15) 



where b m (q) = b- m (—q)*. For simplicity, we impose pe- 
riodic boundary conditions, so that q is an integer mul- 
tiple of 2tt/L. Since the total number of atoms in the 
system is unchanged by the perturbation, 6o(0) is re- 
lated to the other coefficients by volume conservation 



M°) = - 4- £ E IM<z)l 2 + o{\ 2 ), (16) 

' >0 m q>0 



and may be eliminated. Other constraints may be 
utilized to account for confinement potentials more gen- 
eral [37J than the hard walls considered in the present 
article, but do not lead to a qualitative change in the 
stability analysis. 



The energy change (per unit length) under such a 
perturbation is found to be 

^JT = A ' £ £ am{q; R ^ T )\b m (q)\ 2 + 0(\ 3 ), (17) 
m q>0 

where the mode stiffness a m (q) is given by 

a m {q; R, T) = (m 2 - 1)^ + 2n{a s R - ls )q 2 

+ Sa m (q;R,T), (18) 

and Sa m is a mesoscopic electron-shell correction. 

Neglecting for the moment the mesoscopic correction 
fi&m(q)j we find that the perturbation can only lead to an 

instability for m = and qRo < (l — js/agRo) ~ 1, 
which is the criterion for the classical Rayleigh instabil- 
ity pp. Note that a s Ro > 7s for all physically meaning- 
ful radii (c.f. Table Any perturbation breaking axial 
symmetry is classically unfavorable, and we will there- 
fore consider only axisymmetric perturbations (m = 0) 
in the rest of this paper. 

Using semiclassical perturbation theory, the electron- 
shell correction to the mode stiffness for axisymmetric 
deformations was found to be independent of q |2III22| , 

s ^T)=^-^)v sh MR,n (19) 

This turns out to be true only in the semiclassical ap- 
proximation: A fully quantum-mechanical stability anal- 
ysis 001 reveals that long wires undergo a Peierls-type 
instability at q = 2kp , where k^ is the Fermi wavevec- 
tor for subband v. However, the semiclassical results 
are found to provide a good approximation as long as 
the temperature is not too low, and the wires are not 
too long 001- The total mode stiffness ao(q = 1/Ro) 
in the semiclassical approximation is shown in Fig. [21 
together with the density of states <?(£f)- The perturba- 
tion wavevector was chosen so that the surface contri- 
bution to a (dashed curve) is nearly zero. Fig. [3 shows 
that near the thresholds to open new conducting chan- 
nels, where the density of states is large, the wire is 
very unstable (cto < 0). However, in between the sub- 
band thresholds, the shell correction stabilizes the wire 
(ao > 0). 

According to Eqs. fT%|) and fH^I . the most unsta- 
ble mode is m = 0, q = 0. The stability of the wire is 
thus determined by the sign of the stability coefficient 
A(R Ol T) = a o {q = 0;R (h T), 



A(R,T) 



2lXO s 

R 



( d 2 



\dR 2 



ld_ 

RdR 



Vshell{R,T). 



(20) 

For A > 0, the wire is stable with respect to all small 
perturbations, while the wire is unstable for A < 0. The 
stability diagram so determined is shown in Fig. [31 Com- 
petition between surface tension and electron-shell ef- 
fects leads to a complex landscape of stable fingers and 
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Fig. 2 Density of states at the Fermi energy g (top) and 
mode-stiffness for axisymmetric deformations a = ao(q) 
(bottom), normalized by the volume V of the wire. The per- 
turbation wavevector is given by qRo — 1, so that the surface 
contribution (dashed curve) to a is nearly zero. 



arches extending up to very high temperatures: cylindri- 
cal wires whose electrical conductance is a magic number 
1, 3, 6, 12, 17, 23,... times the conductance quantum are 
predicted to be stable with respect to small perturba- 
tions up to temperatures well above the bulk melting 
temperature Tm ~ O.OlTp. This finding suggests that 
metal nanowires are remarkably robust structures. In- 
deed, the principal stable zones shown in Fig. [3] were 
found 26 to persist up to bias voltages eV > O.lep, 
implying that these wires can support current densities 
greater than 10 10 A/cm 2 , which would vaporize a macro- 
scopic wire. (See Sec.|S]for a discussion of nanowire life- 
time as a function of temperature). In Fig. El the values 
|3"2"] a s = £Fkp/80ir and 7 S = AsFkr/iSir 2 , appropriate 
for alkali metals [22], were used. For larger values of a s 
(e.g., for noble metals), the maximum temperatures (in 
units of Tp) of the stable fingers are reduced somewhat, 
but the stability diagram is qualitatively the same. 

The fact, illustrated in Fig. [3| that electron-shell ef- 
fects can overcome the surface-tension driven instability 
of a cylinder is rather remarkable. The surface contribu- 
tion to J? e is O(kpRo), while the electron-shell contribu- 
tion is O(kpRo)^ 1 . For a typical radius kpRo = 10, 
the shell-correction to the energy is thus two orders 
of magnitude smaller than the surface energy! Stability 
is not determined by the energy directly, however, but 
rather by the convexity (or lack thereof) of the energy 
functional, which involves the second derivative with re- 
spect to i?o [cf. Eq. I|2U|)]. Because V s heii is a rapidly 
oscillating function of Rq, its second derivative actually 
has the same characteristic size as the surface contribu- 
tion to the stability coefficient [first term on the r.h.s. of 
Eq. (H3)]. 

Cylinders are special in this respect, because the term 
0(A) in Eq. (|17|) vanishes exactly by symmetry. For a 




Fig. 3 Stability diagram for cylindrical metal nanowires. 
Dark areas indicate stability with respect to small pertur- 
bations, A(Ro,T) > 0. The quantized conductance values of 
some of the stable wires are indicated. 



more general shape (such as a wire with an elliptical 
cross section [2E|) to be stable, the first variations of 
the surface energy and electron-shell energy, which do 
not have the same characteristic size, must cancel. This 
is only possible for small kpRo and/or for small devi- 
ations from cylindrical symmetry. Thus cylinders rep- 
resent about 75% of the experimentally observable 
(most stable) wires, while structures of lesser symmetry 
represent only about 25% of the total. 

Further insight into the stability criterion A(R , T) > 
is provided by the identity 



d^ cy \(Ro,T) 



dR 



Va 

2irR 



A(R ,T), 



(21) 



where ^ cy i(i? ,T) is given by Eq. @ with R(z) = R . 
The wire can lower its free energy via phase separa- 
tion into thicker and thinner segments if and only if 
d^i cy \/dRo < 0. A < therefore corresponds to an inho- 
mogeneous phase [221, while A > corresponds to a ho- 
mogeneous phase. This is confirmed by dynamical simu- 
lation of weakly perturbed stable and unstable cylinders, 
the later evolving into an inhomogeneous wire |33| . 



4 Evolution toward Equilibrium 

In this section, we use the diffusion equation JSJ to study 
the equilibrium shapes of metal nanowires, as well as the 
approach to equilibrium. Figure 0] shows three stages of 
the typical evolution |33U39| of an initially random (a) 
nanowire: After a relatively short time (b), the short- 
wavelength surface roughness is smoothed out, leaving 
a few cylindrical segments, connected by kinks; Eventu- 
ally, all kinks propagate outward and coallesce, yielding 
an equilibrium shape (c) consisting of a cylindrical wire 
suspended between two thicker contacts. 

Several such simulations starting from various ini- 
tial configurations, with conductance ranging from 1 to 
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z 

Fig. 4 Equilibration of an initially random nanowire: (a) 
initial shape; (b) r = 2 x 10 4 ; (c) r = 3 x 10 7 , equilibrium 
structure with G — 12 Go. 

200 G , and lengths 200 < k F L sC 600, all evolved to 
equilibrium structures consisting of one of the stable 
cylinders found in Sec.|31 connecting two quasi-spherical 
contacts (see Fig. EI a)). The shape of the contact is ac- 
tually a close approximation to a Delaunay unduloid of 
revolution j38| , which is a surface of constant mean cur- 
vature, and is an unstable steady state of diffusion equa- 
tion JEJ without the shell-effect term. This is illustrated 
in Fig. comparing the equilibrium wires, rescaled (b) 
by their maximum radius R m ax to a series of unduloids 
(c) of various mean curvature. The curvature of the un- 
duloid is determined solely by the ratio of the radius of 
the cylindrical part to R ma x > and not by the conductance 
of the wire, or its length. In our case, the deep minima of 
the electron-shell potential, Fig. ^ P m the unduloid at 
its connection with the cylindrical part, thus stabilizing 
it. In fact, if one switches off the electron-shell potential 
in the simulations, the equilibrated wires break apart, as 
expected from the Rayleigh instability. The breaking is 
found to happen first at the junction between the cylin- 
der and the lead, suggesting that it is the weak point of 
the equilibrium structure. 

This suggests that the natural evolution of a 
nanowire, at a temperature sufficient for surface atoms 
to diffuse, is to form a cylinder, thus providing an expla- 
nation of the observation of long, almost perfect cylin- 
drical Au nanowires in transmission electron microscope 
(TEM) experiments [8"lll2lll3lll8| . The same type of sim- 
ulation can be used to understand the thinning process 
observed in TEM experiments |19j , where the wire diam- 
eter is seen to decrease step by step through the propa- 
gation of kinks along the wire. 

5 Lifetimes of Metastable Cylinders 

The equilibrium nanowire structures determined in the 
preceding sections are stable with respect to small per- 
turbations, and represent local minima of the free energy 
functional @. However, large perturbations induced 




Fig. 5 (a) Radius R(z) of the equilibrium shapes for 14 sim- 
ulations starting from random initial shapes; the equilibrium 
shapes being symmetric (although the initial shapes are not), 
we only show R(z) for z £ [0, L/2]. (b) Same shapes rescaled 
by their maximum radius Rmax- (c) Series of Delaunay un- 
duloids of various curvature. 

e.g. by thermal fluctuations can drive the nanowire 
out of such a minimum, leading to a finite lifetime of 
these metastable structures. In this section, we use the 
stochastic model derived in Sec. 12. 3. 21 to study this 
process. 

The statistical properties of the stochastically evolv- 
ing field <fi, Eq. fill)) , are described by equilibrium statis- 
tical mechanics. At nonzero temperature, thermal fluctu- 
ations can induce transitions between stable states (i.e., 
local minima) of the potential V ((/)), Eq. (|13fl . Such tran- 
sitions occur via nucleation of a "droplet" of one sta- 
ble configuration in the background of the other, sub- 
sequently quickly spreading to fill the entire spatial do- 
main. When the noise is weak, i.e., at low temperatures 
(compared to the barrier height) most fluctuations will 
not succeed in nucleating a new phase; it is far more 
likely for a small droplet to shrink and vanish. 

A transition state must go "uphill" in energy from 
each stable field configuration. Because of exponential 
suppression of fluctuations as their energy increases, 
there is at low temperature a preferred transition con- 
figuration (saddle) that lies between adjacent minima. 
These are the nucleation pathways. By time-reversal in- 
variance, they are time-reversed zero-noise "downhill" 
trajectories 021- At low temperatures, the expected 
waiting time of the order parameter in a basin of at- 
traction is an exponential random variable, as is typical 
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Fig. 6 Escape barrier AE as a function of the wire length 
L. Here L c is the critical length at which the transition state 
bifurcates. 

of slow rate processes. The activation rate is given in the 
T — > limit by the Kramers formula 

r ~ r exp(-AE/T). (22) 

Here the activation barrier AE is the energy of the tran- 
sition state minus that of the stable state, and To is the 
rate prefactor. The quantities AE and Zo depend on 
the details of the potential, on the length L, and on the 
choice of boundary conditions at the endpoints z = 
and z = L. Based on the equilibrium structures found 
in Sec. 01 we employ Neumann boundary conditions, 
= d z <j){z, t)\ z =a,L- These boundary conditions force nu- 
cleation to begin, preferentially, at the endpoints, con- 
sistent with experimental observations |19j . 

Equation (|14|l with the potential (|13|1 can not in gen- 
eral be solved analytically, but most minima of the po- 
tential V{(f>) can be locally approximated by a cubic po- 
tential 

V ( ± \ ( / } )=-ah + ^l, (23) 

where <f>± = y/a/0 T 4> (<x,0 > 0). The potential 
{V^) biases fluctuations toward smaller (larger) radii. 

Fig. ||J] shows the escape barrier AE as a function of 
the wire length |34|: Below a critical length L c , the tran- 
sition state is a spatially constant field configuration, and 
the escape barrier grows linearly with the wire length L. 
However, at L — L c it bifurcates into a spatially varying 
instanton configuration with characteristic size ~ L Cl so 
that AE becomes length- independent for L ^> L c . 

Our continuum dynamical model thus predicts that 
the lifetime r of a metastable cylindrical nanowire of 
length greater than the critical length L c saturates with 
an escape barrier given by AE^ — liirix—too AE. In 
terms of the physical parameters defining the cubic po- 
tential (GUI, the critical length L c = -^k 1 / 2 / (af3) 1/4 

and AEoo = ^ K 1 / 2 a 5 ^/P 3 / i . The lifetimes t = 1/T 
for several cylindrical sodium nanowires, calculated us- 
ing the best cubic-polynomial fits to the potential 113L 
are tabulated in Table [21 Note that for a wire with 



Table 2 The lifetime r (in seconds) for various cylindrical 
sodium nanowires at temperatures from 75K to 125K. Here 
G is the electrical conductance of the wire, L c is the critical 
length above which the lifetime may be approximated by 
r ~ v^ 1 exp(AE oc /T), and AEoo is the activation energy 
for an infinitely long wire. Note that wires shorter than L c 
are predicted to have shorter lifetimes. 



G 

[Go] 


L c 

[A] 


AE X 
[meV] 


r [s] 


75 K 


100 K 


125 K 


3 


2.8 


250 


4 x 10 5 


2 


5 x 10~ 3 


6 


4.3 


200 


7 


3 x 10~ 3 


3 x 10 -5 


17 


5.0 


260 


7 x 10 5 


3 


8 x 10~ 3 


23 


6.1 


230 


2 x 10 3 


0.2 


9 x 10" 4 


42 


7.2 


250 


2 x 10 5 


1 


10" 3 


51 


6.8 


190 


1 


8 x 10~ 4 


10" 4 


67 


18.8 


180 


0.6 


5 x 10" 4 


7 x 10 -6 


96 


11.4 


250 


10 5 


0.8 


3 x 10 -3 



G/Gq > 1, the lifetime t may not be the typical time 
before the wire breaks, but rather a switching time be- 
tween two different metastable wires with different con- 
ductance values. 

An important prediction given in Table [3 is that the 
lifetimes of the most stable nanowires, while they do ex- 
hibit significant variations from one conductance plateau 
to another, do not vary systematically as a function of 
radius; the activation barriers in Table [21 vary by only 
about 30% from one plateau to another, and the wire 
with a conductance of 96Go has essentially the same life- 
time as that with a conductance of 3Go- In this sense, the 
activation barrier is found to be universal: in any con- 
ductance interval, there are very short-lived wires (not 
shown in Table [2J with very small activation barriers, 
while the longest-lived wires have activation barriers of 
a universal size 

AEoo ~ 0.6 ^ , (24) 
\m e J 

depending only on the surface tension of the material. 
Here m e is the conduction-band effective mass, which is 
comparable to the free-electron rest mass. A compari- 
son of the lifetimes of sodium and gold nanowires |3"4*] 
indicates that gold nanowires are much more stable, as 
expected from the larger value of the surface tension 
cr s (Au) = 5.9cr s (Na). This is consistent with the obser- 
vation that gold nanowires in particular, and noble metal 
nanowires in general, are much more stable than alkali 
metal nanowires. 

The fact that the typical activation energy l|24() is 
independent of Ro may be understood as follows: The 
instanton is a stationary state of Eq. I|12|l : as such, 
the Virial theorem implies that the bending energy 
i^{dz<f>) 2 ) is proportional to {V{4>)). Since k ~ <J s Ro 
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and V ~ l/-Ro, this implies that the characteristic size 
of the instanton L c ~ ^/o^Rq and AEqo ~ \fo~s- 

The lifetimes tabulated for sodium nanowires in Ta- 
ble|21exhibit a rapid decrease in the temperature interval 
between 75K and f 25K. This behavior can explain the 
observed temperature dependence of conductance his- 
tograms for sodium nanowires |11I14I16| . which show 
clear peaks at conductances near the predicted values 
at temperatures below fOOK, but were not reported at 
higher temperatures. 

6 Conclusions 

The NFEM is the simplest possible model of metal 
nanostructures. Nonetheless, it is a remarkably rich 
model, which provides a unified description of quantum 
transport, stability, and structural dynamics of simple 
metal nanowires. ft is hoped that the generic 'proper- 
ties of metal nanostructures elucidated by the NFEM 
can guide the exploration of more elaborate, material- 
specific models, in the same way that the free-electron 
model provides an important theoretical reference point 
from which to understand the complex properties of real 
bulk metals. 
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